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ABSTRACT 


An aerodynamic computer code, capable of predicting unsteady 

and C values for an airfoil undergoing dynamic stall, is used to 
ra 

predict the amplitudes and frequencies of a wing undergoing torsional 
stall flutter. The code, developed at United Technologies Research 
Corporation (UTRC) , is an empirical prediction method designed to yield 
unsteady values of normal force and moment, given the airfoil's static 
coefficient characteristics and the unsteady aerodynamic values, a, A 
and B. In this study, the aforementioned method is Incorporated into a 
dynamics program which, when coupled with these unsteady forces and 
moments, is able to calculate the airfoil's aeroelastic response. Here, 
this aeroelastic program is applied to a three-dimensional, NACA 0012 
wing. The predicted response of this wing is then compared to the 
response of the wing as recorded in a wind tunnel experiment. 

In this experiment, conducted in the PSU 4' x 5 subsonic wind 
tunnel, the wing's elastic axis, torsional spring constant and initial 
angle of attack are varied, and the oscillation amplitudes and 
frequencies of the wing, while undergoing torsional stall flutter, are 
recorded ( 1 .OxlO^<Re^2 . 5x10^ ) . These experimental values show only fair 
comparisons with the predicted responses. Predictions tend to be good 
at low velocities and rather poor at higher velocities. 

These less-than-accurate comparisons are believed to be due to 
differences between ambient, aerodynamic conditions present during the 
collection of the UTRC program's data base and those present during the 
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aforeraenCloned experiment. The major differences are between Reynold s 
Numbers, Mach Numbers and three-dimensional flow effects. It is 
believed tliat eliminating these differences would greatly Increase the 
aeroelastlc program's prediction accuracy. Thus, this program can be 
expected to yield fair predictions for any given set of ambient, 
aerodynamic conditions and probably much better predictions for ambient 
conditions which closely approximate those of the UTRC program s data 
base . 
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CHAPTER I 


INTRODUCTION 


The purpose of this study Is to determine the ability of a 
computer code, capable of calculating unsteady normal force and moment 
coefficients on an airfoil undergoing dynamic stall, to calculate the 
aeroelastic response of a wing undergoing torsional stall flutter. 
Airfoil stall flutter is a periodic oscillation of the airfoil, caused 
by negative aerodynamic pitch damping, which may occur when an airfoil 
experiences incidence angles which vary above and below its static stall 
angle. This phenomenon, observed as long ago as World War I on 
monoplanes, is of great concern in the helicopter and turbomachinery 
industries. It often places operational limits on helicopters and 
turbines by decreasing the efficiency of rotor and compressor blades 
and, in extreme cases, by causing structural failure of these 
components. Recently, stall flutter related problems have been 
encountered on the NASA Advanced Turboprop, a thin, swept propeller 
capable of cruising efficiently at high subsonic Mach Numbers (26). The 
presence of such problems has created the need to predict stall flutter 
onset and subsequent motion of two-dimensional airfoils and three- 
dimensional wings of varied and complex geometries. 

Successful methods of stall flutter prediction have, to this 
date, been restricted to two-dimensional airfoils, and, in particular, 
to the NACA 0012 airfoil or other airfoils of similar geometry. This 
restriction is due mainly to the highly non-linear normal force and 
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pitching moment coefficient behavior characteristic of an airfoil 
experiencing dynamic stall. This characteristic does not allow for 
simplifying linearizing assumptions in numerical model formulations and 
is further complicated by the fact that a quasi-static approach is 
inaccurate, comparable to classical, unsteady aerodynamics. 

A number of dynamic stall prediction methods have been devised 
and, as mentioned earlier, have met with limited success. These methods 
fall into two broad categories, theoretical and empirical. To date, 
empirical methods have given better results than theoretical ones and 
will continue to do so until computer capabilities and numerical methods 
improve enough to allow theoretical predictions to become more reliable. 

The program tested in this study is an empirical one. It was 
incorporated into an aeroelastic code and was then used to calculate 
dynamic, steady-state amplitude and frequency values for a NACA 0012 
wing undergoing torsional stall flutter. These predictions were then 
compared to experimental values obtained from wind tunnel tests done 
with this same wing. The ability of this program to actually predict 
such an aeroelastic response was then determined. 

The main purpose of this paper, then, is to give the reader a 
detailed description and analysis of the outcome of this comparison. In 
addition, the reader will find a description of both the experiment and 
the aforementioned program in Chapters IV and III, respectively, and a 
detailed, but concise, overview of stall flutter and dynamic stall in 


Chapter II. 


CHAPTER II 


AIRFOIL TORSIONAL STALL FLUTTER AND DYNAMIC STALL 
Stall Flutter Characteristics 


Airfoil stall flutter is a predominantly torsional, dynamic, 
aeroelastlc instability caused by negative aerodynamic pitch damping 
which can occur as the airfoil passes in and out of its static stall 
region (10). It differs from classical flutter in that classical 
flutter is dependent upon either elastic and/or aerodynamic mode 
coupling or upon a linear phase lag between the airfoil displacement and 
aerodynamic reaction. It differs from buffet in that the presence of 
buffet forces is solely dependent upon the presence of stall over part 
or all of the airfoil and not upon any subsequent airfoil motion (2). 

In contrast, the negative aerodynamic damping which causes stall flutter 
is a result of the non-linear, unsteady, aerodynamic reaction to airfoil 
motion . 

One of the features peculiar to an airfoil undergoing purely 
torsional stall flutter* is that it oscillates sinusoidally in a limit 
cycle, that is, it oscillates according to the equation 

a = + Aasinwt 


For the remainder of this paper, the term "stall flutter" refers to 
"torsional stall flutter" only, unless otherwise stated. 
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where a-angle of attack 

cx^omean angle of oscillation 
Aa>osclllation amplitude 
u)>osc Illation frequency. 

This happens when the airfoil reaches a dynamic equillhrlum 
where the negative aerodynamic damping from the alrstream equals the 
positive structural damping of the airfoil. Thus, stall flutter can 
occur at virtually any free-stream velocity, V^, where this condition is 
met. Again, this Is In contradistinction to classical flutter where, 
when the flutter velocity Is reached, the airfoil reacts as though being 
forced at a resonant frequency. 

The cause of this limit cycle motion Is Illustrated In Figures 1 
and 2. Notice, In Figure 1, that classical, unsteady, potential and 
values exist when the airfoil Is oscillated through angles of attack 
below Its static stall angle, cigg, and that motion along the 
hysteresis curve as time Increases Is counterclockwise. The exact shape 
of these hysteresis curves Is determined by variables such as reduced 
frequency, k, and oscillation amplitude, Aa (2). The same holds true 
for angle of attack variations which are strictly above the static stall 
angle. However, Figure 2 Illustrates that as the airfoil crosses the 
stall angle, the hysteresis loop begins to distort even to the point of 
crossing over Itself. When this happens, an area Is enclosed which Is 
surrounded by a "clockwise" vs ci hysterep-is curve. This clockwise 
area minus the counterclockwise area equals the net work being done on 
the airfoil by the alrstream and thus represents negative aerodynamic 









damping (6). The magnitude cf this net aerodynamic damping Is Important 
In determining the air full '‘s dynamic response: the larger Its value the 
larger the amplitude and/or the frequency of oscillation. 

It Is possible for any particular airfoil to produce an Infinite 
number of hysteresis loop configurations. To understand why, one must 
have e clear understanding of the mechanism of dynamic stall and the 
factors that most affect Its characteristics. This Is the subject of 
the following two sections. 

Dynamic Stall Characteristics * 

In this section, the major aerodynamic features of a MCA 0012 
airfoil undergoing dynamic stall are first discussed In general and are 
then Illustrated more clearly and quantitatively In an example. All of 
these features are present In varying degrees of relative Importance for 
almost all airfoil geometries and unsteady stalling conditions* Major 
exceptions are noted In this chapter's final section. 

General Characteristics 

C ^ and C^ phase lag (a>0 , This phase lag Is much like 

that which occurs In unsteady, classical potential airfoil theory and 
has been successfully modeled that way (12). 


^Information for the first subsection below Is taken from references 4, 
14-17, 19, 22-24, 27, 31 and 34, collectively, while that for the 
second subsection Is taken solely from references 4, 24 and 30. 
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Stall overshoot . An airfoil undergoing a positive pitch 
velocity will actually "overshoot" its static stall angle without 
stalling* That is, the boundary layer next to the airfoil surface will 

remain attached, and C_ and C will remain virtually unchanged above 

a “h. 

the stall angle. 


Reattachment undershoot . This is the opposite of the stall 
overshoot effect. After the airfoil has gone into full dynamic stall, 
it often "undershoots" its static stall angle before boundary layer 
reattachment begins and before and return to unstalled values. 

Flow reversal and separation . Flow reversal and boundary layer 
separation are not coincident on a dynamically stalling airfoil. In 
fact, a very large portion of an airfoil undergoing dynamic stall may 
experience a thin layer of flow reversal next to its surface before the 
boundary layer actually separates and causes gross changes in and 
values. Figure 3 illustrates tht llfferences in boundary layer profiles 
between static and dynamic separation. 

Vortex shedding . One of the most pronounced events occurlng at, 
or just after, complete boundary layer separation is the shedding of a 
vortex which travels from the leading edge to the trailing edge of the 
airfoil. This vortex, which is usually shed at or very near the maximum 
angle of attack, is the cause for the characteristic extreme values of 

^n-max ^m-mln* 



a. Steady separation 



FIGURE 3: Boundary layer profile^, for static and dynamic airfoil stall 
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Separate stall events . "stalls" just before C„ does, due to 
the effect of the shed vortex mentioned above. 

Force and moment hysteresis . and C„, when graphed as a 
function of a, form distorted hysteresis loops (see previous section). 

Typical Cycle 

The following paragraphs describe a NACA 0012 airfoil undergoing 
sinusoidal pitch oscillations about its 1/4 chord. Pertinent data is as 
follows: a -14.0°, a-15°+10Oslna)t, k-0.15, Re-2. 5x10^ and chord-4 ft. 

static stall ang] - . exceeded . The airfoil, pitching up, behaves 
as though it is still below its static stall angle. No flow reversal is 

present . 

Flow reversal begins (Figures 4a and 4b). When the angle of 
attack is equal to 19 or 20 degrees, a thin layer of reversed flow 
appears at the trailing edge and moves up toward the leading edge. The 
boundary layer remains attached but is thickened and disturbed. No 
noticeable change in or occurs until the flow reversal reaches 

x/c-0.50. 

Vortex formation . This flow reversal continues toward the 
leading edge until it reaches x/c-0.30 (ct-23.4°), at which time the 
boundary layer abruptly separates all the way to the leading edge. 
Immediately, as shown in Figure 4c, a vortex begins to form at x/c-0.1, 
indicated by Increased suction at this point on the airfoil, and moves 
downstream at approximately 35 to 40 percent of the freestream speed 
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(Figure 4d). As this vortex travels, the magnitude of the reversed flow 
near the airfoil surface Increases. 

C and C_ changes. This change in velocity, and therefore 

— n m — — * — 

pressure distribution, over the leeward side of the airfoil causes gross 
changes in force and moment values. For instance, the changing pressure 
distribution causes to "stall," and vortex suction causes to 
increase far beyond its unstalled value. reached when the 

vortex is at x/c-0.50, as shown in Figure 4e, after which begins to 
stall, reached just before the vortex passes over the 

airfoil's trailing edge. 

Reattachment . After the passing of the vortex, the airfoil is 
fully stalled (Figure 4f), l.e., the boundary layer is fully separated 
and Cjj and are in their stalled non-linear range. Clockwise C^j 
hysteresis is greatest at this point. Reattachment begins on the 
downstroke when a»16.5° and progresses rearward from the leading eoge at 
approximately 25 to 35 percent of the free-stream velocity (Figure 4g). 
Full reattachment does not occur until a»7.0°, and unstalled and Cj^ 
values do not occur until a-6.0'' on tha upstroke. The cycle is then 


repeated. 
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Factors Affecting Dynamic Stall 

The above description is typical of an airfoil undergoing, what 
is known as, deep dynamic stall. Deep dynamic stall is roughly defined 
as dynamic stall in which unsteady force and moment values deviate 
greatly from their static counterparts. Very generally, it occurs when 
the maximum angle of incidence is more than Just a few degrees above 
and/or when the reduced frequency is at, or above, a value of 
approximately 0.10. All airfoils undergoing deep dynamic stall exhibit, 
qualitatively, the same characteristics, w!ilch are nearly Independent of 
airfoil shape, Reynolds Number, and type of motion. This consists of 
large fluctuations in and and a well defined and powerful, shed 
vortex. The other "type" of dynamic stall, known as light dynamic 
stall, is characterized by a relatively small deviation of unsteady 
force and moment values from their static counterparts but rather large 
hysteresis and net negative aerodynamic damping. This negative damping 
is usually greater than that which occurs in deep dynamic stall. Vortex 
activity and influence are also considerably less and static 
characteristics are generally much more dependent upon all of the 
factors listed below. 

The purpose of this section is to briefly describe the major 
factors affecting the occurrence and relative importance of each of the 
aforementioned dynamic stall characteristics as they pertain to the NACA 
0012 and its modifications*. 

*For details concerning more advanced airfoil sections see reference 30. 
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All foil Geometry 


in the light dynamic stall region, airfoil geometry haa a 
greater effect on dynamic stall events than any of the other factors 
discussed in this section. Ceometrlcal specifications, such as leading 
edge radius and camber, are usually responsible for the airfoil's 
characteristic static boundary layer separation (eg.. le»)lng edge, 
trailing edge, thin airfoil etc). The nature of the static separation 
greatly affects the nature of the dynamic separation and therefore, 
affects dynamic stall characteristics. Houever, once an airfoil 
penetrates deeply Into the stall region, leading edge separation 
predominates, regardless of airfoil geometry. 

.educed leading edae radius . As Figure 5 demonstrates, there 
are almost no changes In dynamic or static stalling characteristics for 
this slight modification of the NACA 0012. The NACA 0012 stalls with a 
rather abrupt trailing edge separation. 

.ham leading edge (Figure 6). In this configuration, both 
static and dynamic stall begin with laminar leading edge boundery layer 
separation. During dynamic stall this cause, flo. reversal to begin 
earlier and more abruptly. In addition, vortex formation and movement 
is coincident with flow reversal Incipience and spreading. In ocher 
words, they begin slmulteneously at the leading edge and travel together 
toward the trailing edge. These differences cause significant changes 

in unsteady and curves. 
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Cambered Leading Edge . Here, both static and dynamic stall 
remain of the trailing edge type but they become much more gradual. The 
shed vortex Is much less powerful, and, as shown In Figure 7, this 
results In lower and values. 

Reduced Frequency 

Next to airfoil geometry, reduced frequency, ujc/2V^, has the 
most significant effect on dynamic stall. Generally, Increases In k 
cause both a delay of stall events, such as flow reversal, boundary 
layer separation and vortex formation, and an Increase In the phase 
angle between them. Specifically, an Increase In k usually causes a 
delay In the occurrence of flow reversal . This delay, in turn, causes 
separation delay and a proportionate delay in vortex shedding. A vortex 
which Is shed later is generally more powerful due to the presence of an 
increased attached circulation before stall. Consequently, C^-^ax 
C , become greater; their values can change by as much as 100 percent 
as k is varied from 0.004 to 0.25, where k^O.004 is generally considered 
to be the quasi-static region. Figure 8 shows how these unsteady values 
vary with k. In addition, a similar delay In boundary layer 
reattachment causes larger hysteresis and, up to a point, greater net 
negative aerodynamic damping. As stated earlier, when conditions for 
deep dynamic stall are reached, net negative aerodynamic damping begins 
to drop off despite large and values. Thus, increasing k 

*It has been proposed that this delay Is due to the filling out of the 
boundary layer near the airfoil's leading edge, resulting from the 
positive pitch velocity (12). 
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will cause greater differences between static and dynamic stall 
behavior. 

Amplitude and Mean Angle 

These two factors are important in that they define the 
incidence angles through which the airfoil sweeps. This, in conjunction 
with reduced frequency, determines whether the airfoil experiences light 
or deep dynamic stall. A large and reduced frequency will generate 

a strong bound vortex, due to stall delay at high incidence angles and a 
well defined and powerful shed vortex. In contrast, an 
greater than will cause either no stall at higher values of k. as 
illustrated in Figure 9, or light stall at lower values of k. Again, 
amplitude and mean angle, along with reduced frequency, take on greater 
importance than airfoil geometry during deep dynamic stall oscillations 
due to the fact that under such conditions the boundary layer tends to 
separate near the leading edge of the airfoil. 

R eynolds Number and Mach Number 

Reynolds Number plays a rather insignificant part in changing 

dynamic stall characteristics, at least between the common test values 
of 1.5x10^ and 3.5x10^. However, it has been found that, at higher 
values of Reynolds Number, flow reversal becomes more abrupt while, at 
lower values, a greater moment hysteresis results due to a delay in flow 
reattachment. This can possibly cause significant increases in net 
negative aerodynamic damping at very low Reynolds Numbers. 
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Mach Numbac generally haa little or no effect on dynamic stall 
• as long as It la kept below a value of approximately 0.2. However, 

I above tills value, but below a value of 0.6, and begin to 

decrease due to premature, shock Induced separation. Above M -0.6. no 
I overshoot of static valuea la seen though hysteresis still may occur. 

In addition. Increased subsonic Mach Numbers, like deep stall 
* conditions, cause unsteady separation to tend toward the laminar leading 

I edge type regardless of airfoil geometry. 

I 

I 
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CHAPTER III 


COMPUTER PREDICTION OF STALL FLUTTER 

A program that can accurately predict unsteady C^^ and Cj^^ values 
for an airfoil undergoing dynamic stall, coupled with a dynamics 
subroutine, should be capable of doing a reliable stall flutter 
analysis. However, the difficulty In predicting these highly non- 
linear, unsteady values limits the accuracy of this type of program. 
Despite this difficulty, a number of programs exist which are capable of 
yielding fairly good results. In the sections below, the reader will 
find a brief description of the existing types of dynamic stall 
prediction methods and a detailed description of an empirical method 
developed at United Technologies Research Corporation (UTRC). In 
addition, the reader will also find a detailed description of how this 
UTRC dynamic stall prediction method was Incorporated Into the stall 
flutter prediction program that was used In this study for comparison 
with experimental data. 


Theoretical Models 

A detailed description of theoretical, dynamic- stall prediction 
programs Is beyond the scope of this paper; however, their basic 
features, as reported In references 25 and 28, are outlined below along 
with a discussion of the difficulties currently being encountered In 
making these programs practical. 

A good numerical modex, short of one that solves the full non- 
linear Navler-Stokes Equations, should Include elements of the 


- 23 - 


following: unsteady boundary layer theory, potential vortex theory and 
viscous-inviscid flow regime interaction. The problem with this type of 
approach is that the program does not "understand" the difference 
between unsteady flow reversal and separation and that it cannot predict 
the type of boundary layer separation that will occur (e«g.» leading 
edge, trailing edge etc.). In addition, it must also be given empirical 
estimates of vortex incipience, strength, velocity and viscous 
dissipation. Add to this inaccuracies still present in complex flow- 
field modelling and predicted results become relatively unreliable. 
Probably the beat way to theoretically predict dynamic stall is to solve 
the full Navier-Stokes Equations. Unfortunately, computers are not yet 
large enough or fast enough to handle this problem. 

Empirical Methods* 

Empirical programs, concerned with predicting dynamic stall, 
relate unsteady airfoil normal force and moment data to their static 
counterparts. This is usually done by using dimensionless parameters, 
such as o)c/V^, db/V„, db^/y„^, coupled with empirical functions designed 
to take unsteady stalling characteristics into account. The data base 
fcr these types of programs is most often obtained by oscillating an 
airfoil sinusoidally about its 1/4 chord over a large range of values of 
reduced frequency, mean angle of attack and amplitude. 


*For a concise overview of these methods see reference 27. Individual 
methods can be found in references 1,3, 9-13, 17 and 20. 
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In general, such an approach has been rather successful In 
reproducing basic dynamic stall characteristics (eg*, lift stall 
overshoot, loop hysteresis etc)* Not only do these programs reproduce 
their own data base well, but they have been shown to be fairly accurate 
when predicting a set of Independent sinusoidal data (28)* This point 
will be discussed In greater detail later In this chapter* In addition, 
they are self-contained; except for static data curves, few. If any, 
other data need to be Input Into the program before execution* 

However, the general applicability of an empirical method to an 
airfoil geometry which 1s dlf fete-, t than that of the airfoil used to 
obtain the data base Is questionable* This could be particularly 
important in the light stall region where airfoil geometry greatly 
affects stall characteristics and where, incidentally, most practical 
cases of stall flutter occur* 


United Technologies Djmamlc Stall 
Prediction Program 

Data Base (5) 

The program chosen for comparison with experimental data. In 
this study. Is an empirical one developed at UTRC in 1975* This program 
uses the parameters a, A and B to calculate unsteady and Cjjj values 
from their steady counterparts* The data base was collected from a NACA 
0012 airfoil, having a span of 3*75 inches and a chord of 5*0 Inches, in 
the United Aircraft Corporation's two-dimensional, high-speed, cascade 
wind tunnel (test section: 4” x 14” high)* The airfoil was oscillated 
through sinusoidal, ramp and, what Is called, helicopter angle of attack 
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variations at a free-atream Mach Number equal to 0.325 and a Reynolds 
Number equal to 9.6x10^. Parametric variations for data collection are 
shown In Table I. After being collected, these unsteady and Cj^ data 
were tabulated as a function of the Instantaneous values of a, A and B. 

Curve-Fitted Equations 

In 1975 these thousands of tabulated data points were placed In 
a more convenient format by curve fitting them, again according to a, A 
and B, using a non-linear, least-squares method (I). This was done by 
"expanding" the static and (vs d) lines into the dynamic ones in 
the following way: 

1. The linear parts of the static curves were matched to the 
linear parts of the dynamic curves through a horizontal and a 
vertical shift of the static stall point, as shown in Figure 10. 
The equations that affect the horizontal shift, Aa, are 


(Aa/a )„ = A(P„,+P„£a,)+B(P_+P„, a,) 
ss N N1 N5 N N2 N6 N 

,2 




^- 40316.428 


and 






SS 


where o»a/a, 
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TABLE I 

TEST PROGRAM OUTLINE 



Cam Description 


Primary 

Drive 

Mean Angle* 



[tern <?nd Angular 

Type of 

Frequency 

of Attack 



No. Amplitude 

Coupler 

(cps) 

(deg) 

—— 1 


3 

4 


10 


11 


noue 


rigid 


0 (steady state) 0,3,7,9,10, 






11,12,13,14, 

16,18,20,22 

1 

sinusoidal. 

± 8 

rigid 

12.5,31,50,75 

100 

3.6,9,11,12, 

14,16,18 ! 

1 

sinusoidal. 

± 8 

rigid 

120 

1 

3 1 

sinusoidal. 

+ 8 

rigid 

50,75 

1 

11 (single 
surface meas.) 

sinusoidal. 

± 3 

rigid 

10,25,50,75 

14 

helicopter. 

± 8 

rigid 

10,16.7, 33 

6,11,14 

forward ramp, 
#8 

rigid 

7.5,10,15,20 

6,11,14 

1 

1 

1 

backward ramp, 
#8 

rigid 

7.5,10,15,20 

6,11,14 1 

1 

1 

sinusoidal. 

#8 

flexible 
(3 springs) 

7.5,10,15 

1 

6,11,16 1 

1 

1 

1 

none 


flexible 
(3 springs) 

0 (free 
vibration) 

6,11,16+a 1 

for max. | 

oscillation | 

1 

none 


flexible 
(3 springs) 

impulsive 
loading, Aas3 

1 

6,11,16 1 

1 

1 


*Note: Each item represents test point taken at all possible combinations 
of frequency and angle of attack. 
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TABLE II 

. COEFFICIENTS for SYNTHES IZ AT ION 
FUNCTION EVALUATION. 

•f 

1 

1 

1 

1 




* 


1 

1 1 1 

^N1 j 

Qni 

I 

^Mi 

1 

1 

1 

. 



. 


, ^ 




• 



1 1 1 

-3.464003-01 1 

1.533717+00 

1.970065+01 1 

-2.322 808+00 

1 

1 2 1 

-1.549076+00 ! 

6.977203+00 

-6.751639+01 1 

-1.322257+00 

1 

1 3 1 

4.306330+01 1 

1.749010+03 

7.2652 69+02 | 

-2.633891+00 

1 

1 4 1 

-5.39752 9+01 I 

1.694829+03 

4.865945+04 | 

-2.180321-01 

1 

1 5 1 

5.7814QZ+00 1 

-1.771899+03 

2.0862 79+04 | 

4.580014+00 

1 

1 6 1 

-3.2 33003+01 1 

-3.291665+04 

6.024672+03 I 

3.125497-01 

1 

1 7 1 

-2.162257+01 1 

2.969051+00 

1.446334+02 | 

-2.82 8806+01 

1 

1 8 1 

1.866347+01 1 

-3.632 448+01 

8.586896+02 j 

-4.396734+00 

1 

1 9 1 

4.198390+01 1 

-2.268578+03 

-7.550329+02 | 

2.565870+02 

1 

1 10 1 

3.295461+02 | 

6.601995+03 

-1.021613+01 1 

-1.204976+01 

1 

111 1 

— 1 

-9.6542 08+03 

2.247664+01 1 

-1.157802+02 

1 

1 12 1 

— 1 

8.533930+04 

— 1 

8.612138+00 

1 

|13 1 

— 1 

-1.492 624+00 

— 1 

•— 

1 

|14 1 

1 

1.163661+01 



1 

•f 


and where the values of P and Q are given In Table II. 
The vertical shift is executed by multiplying the linear 
coefficient slope by the appropriate value above. 


2. Residual differences in both the linear (ACj) and non-linear 
(AC 2 ) range are "set to zero" with the following equations 

AC^=A[Qni+Qn 3A +C^(Qj^^+QjjgA ^%5‘‘^N11%^^ 

B[Qn2'^N4^ +<^(Qjj8+Qji10^ ^^N6‘^N12 ^ ’ 




3> These values are then combined Into one equation to 
determine the final unsteady value, 

u N 

u M 

where a^-Unear coefficient curve slope with respect to a. 

4. Scaling for non-NACA 0012 airfoils and for Mach Numbers 
other than 0.325 Is rather straightforward. Different 
geometries are accommodated by reading the static data (0^^ and 
Cjjj vs a) for that particular airfoil Into the program before 
execution begins. Mach Number effects are taken Into account by 
using static data obtained at that particular Mach Number and by 
multiplying AC^ and by the Prandtl-Glauert Correction 
Factor, 6, for compressible flow and by the constant 0.9457. 
Thus, the final equations become, with corrections, 

n Q4S7 

" ** "W 

5, <o,A.B.M)=.C„(c.A.B.M) .a (M)Ac^f 

“ M W 
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As seen by Figures 11a through 11c, these equations do well at 
^®P^- 0 ^^clng their own data base* Exact data points are usually not 
predicted but, again, basic dynamic stall characteristics are predicted* 
There are definite bounds on a, A and B, defined by the data 
base limits, within which this program remains valid* To get a picture 
of how these limits are defined , Imagine values of a , A and B graphed In 
three-dimensional space* They form an elliptic cylinder (A and B are 
semi-major and semi-minor axes, respectively) with irregular endplates, 
a* Now, Instead of limiting unsteady values of A and B to those 

by this ellipse, they are. Instead, defined aA follows: 

A,B n£i 

(A,B) = 

i(A,B) n>i 

This will automatically transfer any point (A,B) outside of the ellipse 
to a point Inside of the ellipse* Values of a are, however, restricted 
to their extrema: if a (a/oigg) becomes greater than either 1*86 or 
1* 839-70 *33B, the program becomes Invalid* This limits a to maxima and 
minima under various unsteady condititions* 
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Independent Evaluation of Empirical Methods 

To date, there has been only one comparison between predictions 
obtained from empirical programs and forced sinusoidal data collected 
from an Independent source. The outcome of this comparison Is reported 
In references 28 and 29. This was done to check the ability of these 
methods to predict data from a source other than their own data base. 
Figures 12a and I2b, taken from reference 28, Illustrate just how 
closely five existing empirical programs predicted values of 
^m-min' which were taken from an Independent set of NACA 0012, dynamic - 
stall data. Results from this comparison showed fairly good agreement 
for all methods tested. Notice that the UTRC program showed 
consistently low predictions of both and . This may be due 

to the relatively high Mach Number at which the UTRC data base was 
collected. Shock-Induced stall may well have caused early separation 
and, consequently, smaller extreme values for and as explained in 
Chapter II. 


Airfoil Aeroelastlc Response 

With a means available for calculating unsteady, non-linear 
normal force and moment values, the next step Is to Incorporate these 
Into a dynamics program to predict the total aeroelastlc response. A 
program with this capability has been run on the PSU IBM 370/3033 
Interactive computer system*. Details as to how this has been dcue are 

^Copies of the program listing are available upon request. 



FIGURE 12: Comparison between empirical program predictions and ^ 
independent data base. (oi=15°^°sinojt, k*0.15, Re-2. 5x10 ) 


UTRC I UTRC 








I 

I 

I 

I 

I 

I 




b. 

FIGURE 12 (cont.)* (ct=15°+10°sincot , k=0.15, Re=2.5xl0^) 
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given below. 

An airfoil is placed in a free-stream velocity with initial 
conditions on the torsional displacement angle, 6, and angular velocity, 
g, specified, as shown in Figure 13. Unsteady coefficient values are 
computed by the UTRC program, assuming initial values for a, A and B, 
and are used to "drive" the airfoil dynamics. The motion of an airfoil 
undergoing single-degree-of-freedom torsional oscillations is described 

by the equation 

le I3+K06 -(Cm+(EA-AC)Cl)xO . spy^^sc , 

where I=-moment of Inertia of the wing about its elastic axis 
^■damping ratio 

M -natural frequency of the wing 
n 

and where other variable definitions are given in Figure 13. 

Specifically, this non-homogenou.- , second-order differential equation is 

solved by the fourth-order Runge Kutta method over a short time 

interval, during which the unsteady coefficients are assumed to remain 

constant (18). This solution yields new "initial conditions" for the 

unsteady parameters, and these are used to calculate updated values of 

C and C . These updated values again drive the dynamics, and the cycle 
n m 

is repeated until dynamic equilibrium is reached. Figure 14 shows the 


flow diagram for this program. 






FIGURE :.4: Flow diagram for aeroelastic program 







CHAPTER IV 


I 

I 
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STALL FLUTTER EXPERIMENT 

The purpose of the experiment described In this chapter was to 
collect dynamic data (amplitude and frequency vs. velocity) for a NACA 
0012 wing section undergoing torsional stall flutter. This experimental 
data was used as a standard with which to compare data predicted by the 
aeroelastlc program described In Chapter III. Results of the comparison 
are given In the next chapter. 


Apparatus 

This experiment was conducted In the PSU 4' x 5' subsonic wind 
tunnel, a closed circuit, atmospheric tunnel capable of a maximum 
velocity of approximately 140 ft/sec. The wing model was a hand-crafted 
rectangular, unswept NACA 0012 made of solid Phllllplne Mohogany, coated 
with about 15 layers of lacquer. It had an aspect ratio of 5.3 and a 
chord length of 4.5 Inches. Circular endplates, 1 foot In diameter and 
made of 1/8" plexiglass, were attached to the wing tips. 

Static data for the wing, which was measured with a pyramidal 
balance placed below the tunnel test section, was taken with the 
endplates, aluminum support bars and accelerometer connected (see next 
section), since these attachments all rotated with the wing during 
flutter. This data, shown in Figures 15 and 16, is uncorrected for 
solid blocking, wake blocking or streamline curvature (32). It was left 
in this "raw" form for two reasons: first, these effects were estimated 
to be rather small, due to the relatively small size of the wing as 
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FIGURE 16: Moment characterist 
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compared to the test section, and second, the dynamic data was taken In 
the same tunnel where the same effects would be present. Drag 
Interference was corrected, to a first order approximation, by taking 
measurements on the static-measurement-support bar without the wing 
attached and subtracting these values from those obtained with the wing 
attached . 

The basic shapes of these static curves are like those of a NACA 

0012 except that the stall angle Is slightly below 11.0 degrees. This 

Is most likely a Reynolds Number effect (2.5x10^). Notice the abrupt 

stall of the lift curve, characteristic of a 0012. Tuft visualization 

verifies that this stall Is of the abrupt, trailing edge type, again 

characteristic of the 0012. Finally, the lift curve slope, 

approximately equal to 0.0914/degree, Is a reasonable value for this 

geometry, but the value of at zero angle of attack of 0.015 should, 

of course, equal zero for a synmetrlcal wing. This discrepancy. In 

addition to the stray points and positive , evident on the moment 

a 

curve. Is probably due to a combination of aluminum support bar and 
accelerometer contributions and Inaccuracy In the airfoil construction*. 

The frame used to hold the wing was made of 1/8” x 1 1/4” flat 
steel bars, welded together In a rectangular-box configuration (see 


A repeat of the static tests verified the data shown In Figures 15 and 
16. However, It should be noted that while general features, such as 
stall angle and general curve shape, were repeated In the re-testlng, 
specific data points were rarely repeated. In addition, the updated 
static characteristics consistently yielded a slope above 
0.1/degree, a value greater than a NACA 0012 airfoil section. This 
should never occur for a three-dimensional wing. The reason for this 
large llfi: curve slope Is unknown. 
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Photograph 1) . It had an overall height of 20.5" a "span" (perpendicular 
to the wind velocity) of 27.0" and a "chord" of 10.0 . A hinged wooden 
bar, connected to the tunnel floor and downstream legs of the frame, 
enabled it to rotate to varying angles of attack. It was "locked out" 
at these different angles in the following way: two threaded rods, 

bolted to the tunnel floor, "pierced" the "eyes" of two eyebolts which 
were attached to the frame's forward legs. These eyebolts were then 
held into position with two opposing nuts mounted on each rod. The wing 
and endplates were attached to the frame, approximately 15.0 above the 
tunnel floor, via two 16.0" aluminum bars and six 1 1/2 wood screws. 

The torsional restoring spring consisted of between 2 and 8 
translational linear springs*. They were 1 1/2" long and were attached 
to the frame and the aluminum bars approximately 7.5 from the wing s 
center of rotation. The elastic axis was varied, then, by positioning 
the wing such that the midpoint of the aluminum support bar was adjacent 
to the desired elastic axis. Finally, an electromagnet was hung from 
the top chord -beam on the frame's far side. It was used to prevent the 
wing from oscillating until the tunnel wind could be brought up to the 
desired speed. 


*See Appendix A for calculation of torsional spring constants. 
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Experimental Procedure and Data Measurement 


When the desired torsional spring constant, Initial angle of 
attack, and elastic axis were set, the electromagnet was turned on, and 
the wind velocity was brought up to speed (test Re: 1.0x10^-2.5x10^). 

The magnet was then switched off and, if flutter occured, the free- 
stream velocity, frequency, peak-to-peak acceleration, and maximum and 
minimum angles of oscillation were recorded. If, when the magnet was 
switched off, flutter did not occur, the wing was tapped until it did 
occur. The maximum angle to which the wing had to be tapped, in order 
to excite self -sustained flutter, was then recorded along with the above 
mentioned data. After this, the tunnel velocity was brought to zero, 
and the procedure was repeated for a different value of V^^,. These steps 


I TABLE III. VARIATION of PARAMETERS 
I for PSU EXPERIMENT. 

I 


x/c 

VELOCITY 

(ft/sec) 

"i/%s 

L . J 

‘"n 

(cyc/sec) 

L . . 

0.250 

28-110 

0.64-1.27 

1.85-3.82 

0.375 

28-70 

0.45-1.36 

3.22-6.0 

0.500 

28-70 

0.64-1.36 

3.30-5.90 


were carried out for the parametric varations shown in Table III. 

Static equillbrian angles were measured using angular markings 
drawn on a clear plastic surface which was taped to the near end of the 
frame. Dynamic flutter amplitudes were measured using both the 
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aforementioned markings and an Endevco 2233 Accelerometer, which was 
attached to the near-end, aluminum wing-support bar. The 
accelerometer's signals were wired through an Endevco Amplifier and a 
high-frequency filter before being displayed on an oscilloscope screen 
from which peak-to-peak amplitude and frequency values were taken. 

Experimental Results 

Limit cycle oscillations about the 1/4 chord elastic axis were 
found to be very nearly sinusoidal. Flutter amplitudes varied from 
approximately 2.0 to 7.0 degrees. In general, as can be seen from 
Figure 17, the amplitudes increased with respect to V^. Values for uj/uJjj 
varied from 1.1 to 1.3 and showed a slight increase with respect to , 
as shown in Figure 18. In addition, k varied from 0.12 for stiff 
springs at low V^, to 0.03 for "loose" springs, at high values large 
enough to yield unsteady stalling for all conditions studied. It should 
be noted here that a similar experiment, done at MIT (8), showed 
amplitudes of the same magnitude but produced no increase in these 
values with V^, for purely torsional stall flutter. However, this 
experiment did yield an increase in w with V^. Figures 19 and 20 are 
graphs of these values. 

Limit cycle oscillations about the 37. 5Z chord elastic ai.is were 
also found to be sinusoidal, but the amplitudes, as expected, were much 
larger, ranging from approximately 4.0 to 15.0 degrees. Again, as 
Figure 21 shows, amplitudes increased with but began to level off at 
higher velocities. However, for this case, stayed equal to 1.0 for 
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FIGURE 17: (cont.). 
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FIGURE 18: Experimental frequency values for PSU wing undergoing 
torsional stall flutter. 


Elastic Axis*0.250c 
Natural Frequency*2 . 1 0 Hz 
Initial Angle of Attack=H.0 Degrees 
Zeta=0.01050 





]pTTnrnpTmn n [i n iMt t TpTrmtfTpnTim|TmmqmmrTTjntTMiii|inii.iii|iiMiiiii|iiiiiniijiiiiimT|- 

0 10 20 30 40 50 60 70 80 90 100 110 120 


VaOCITY (ft/sec) 


b. 

FIGURE 18: (cont.)» 
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FIGURE 18: (cont.). 
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FIGURE 18: (cont.). 
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FIGURE 19: Experiir-ancal amplitude values for MIT airfoil undergoing torsional 
stall flutter. 
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FIGURE 20: Experimental frequency values for MIT airfoil undergoing 
torsional stall flutter. 
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FIGURE 21. Experimental amplitude values for PSU wing undergoing torsional 
stall flutter. 
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FIGURE 21; (cont.). 




Elastic Axis=0.375c 

Natural Frequency^S . 83 Hz 

Initial Angle of Attack=10.5 Degrees 

Zeta=0.00888 



0 10 20 30 40 50 60 70 80 


VaOCITY (ft/sec) 


c. 


FIGURE 21: (cont.). 



62 


Elastic Axis=0. 375c 

I Natural Fr equency=6. 00 Hz 

Initial Angle of Attack=9.0 Dagreos 

I 2eta=0.00493 

1 20j 



0 18 20 30 « S8 60 70 80 30 100 

VaOCITY Cft/sec) 


j 

( 


figure 21: (cont.). 
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all trials and conditions (Figure 22). Again, reduced frequencies were 
high enough to produce unsteady conditions; values of k ranged from 0.04 

to 0.23. 

Results for oscillations about the 50.0% chord elastic axis are 
virtually the same as those obtained from oscillations about the 37.5% 
chord elastic axis, except that amplitudes vary from approximately 10.0 
to 20.0 degrees, as shown in Figures 23 and 24. 
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FIGURE 22: Experimental frequency values for PSU wing undergoing torsional 
stall flutter. 
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FIGURE 22: (cont.). 
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FIGURE 22: (cont.) • 
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FIGURE 22; (cont.). 
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FIGURE 23; Experimental amplitude values for PSU wing undergoing torsional 
stall flutter. 
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Elastic Axis=0.ffl0c 
Natural Frequency=3 . 83 Hz 
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FIGURE 23: (cont.). 
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FIGURE 23: (cont.). 
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FIGURE 23: (cont.). 
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figure 24= Experimental frequency values for PSU wing undergoing torsional 
Stall flutter. 
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FIGURE 24: (cont.)- 
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FIGURE 24: (cont.). 
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FIGURE 24: (cont.). 


CHAPTER V 


COMPARISON BETWEEN PREDICTED AND EXPERIMENTAL 
STALL FLUTTER RESPONSE 

The comparison between predicted and experimental results is 
only fair. For the case where the elastic axis was positioned at the 
1/4 chord, predicted values of amplitudes are consistently lower than 
those taken from the experiment. Figures 25a through 25d show typical 
cases. Note that comparisons are fairly accurate at lower velocities, 
but that they get progressively worse at higher velocities. This is 
because computed values show little increase with , in contrast to the 
similar but much greater increase shown by experiment£il data. Relative 
values of frequency show exactly opposite trends: predicted values of cd 

overestimate experimental values, and they Increase with V^, while 
experimental values do not (Figure 26). Comparison with MIT data, where 
MIT's NACa 0012 airfoil data were substituted into the UTRC program, 
shows an underprediction of both amplitude and frequency, but the slopes 
of both are accurately predicted (Figure 27). 

Predicted and experimental data taken from the 37.5% chord 
elastic axis show much better agreement than those taken from the 1/4 
chord elastic axis. Unfortunately, the program often went out of its 
bounds due to the high angles of attack reached during the larger 
amplitude oscillations. When it did, predicted values diverged yielding 
very poor comparisons. However, when it stayed in bounds, predicted 
values yielded the following results: amplitude errors are of the same 

order of magnitude as those of the 1/4 chord elastic axis, but, in 
general, the slope of the amplitude lines are predicted fairly well 
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FIGURE 25: Comparison between PSU experimental amplitude values and 
values predicted with the UTRC program. 
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FIGURE 25; (cont.). 
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FIGURE 26: Comparison between PSU experimental frequency values and values 
predicted vith the UTRC program. 
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FIGURE 26 : (cont.). 
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FIGURE 26: (cont.). 
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figure 27: Comparison between MIT experimental amplitude values and 
values predicted with the UTRC program. 
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FIGURE 27; (cont.). 



(Figure 28). However, in contrast to 1/4 chord results, program 
amplitude predictions overestimate experimental values. Finally, 
frequency values are predicted almost exactly, as shown in Figure 29. 

Experimental oscillation amplitudes about the 50. OZ chord 
elastic axis are so large that predictions calculated at only the lowest 
velocities, yielding the smallest amplitudes, are in bounds. 

Prediction of stall flutter onset was good in some cases and not 
so good in others. In most cases, predicted stall flutter did not begin 
on its own unless the experimental flutter did (le., unless experimental 
stall flutter began without having to tap the wing), and vice-versa. 

Most of the discrepancies were encountered when predicting the angle to 
which the wing had to be tapped. For Instance, a fairly good prediction 


TABLE IV. FLUTTER ONSET PREDICTION 


VELOCITY 

(ft/sec) 


^1-e 

(deg) 


WING RESPONSE 
(experimental) 


FLUTTER 


FLUTTER 


FLUTTER 


^1-p 

(deg) 


WING RESPONSE | 
(predicted) | 

FLUTTER I 


FLUTTER 


FLUTTER 


is shown in Table IV. Here, it is found that, when uas set equal 
to self-sustained stall flutter occurred, whereas, without this 

no flutter occurred. A poor prediction is illustrated in 
Table V. Gross Inconsistencies exist here. Notice, that at 46.3 ft/ sec 
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FIGURE 28; Comparison between PSU experimental amplitude values and values 
predicted with the UTRC program. 
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FIGURE 28: (cont.)« 
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figure 28 : (cont.). 
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FIGURE 29: Comparison between PSU experimental frequency values and values 
predicted with the UTRC program. 
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FIGURE 29; (cont.). 
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FIGURE 29: (cont.). 
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TABLE V. FLUTTER ONSET PREDICTION 


+ 


VELOCITY 

(ft/aec) 


(deg) 


WING RESPONSE 
(experimental) 


®i-P 
I (deg) 


WING RESPONSE 
(predicted) 


4 - +_ 1 + 


46.3 

1 

00 1 
1 

NO FLUTTER 

1 

1 

10.0 1 
1 

FLUTTER 

55.1 

1 

1 

1 

5.5 1 
1 

FLUTTER 

1 

1 

1 

7.0 1 
1 

NO FLUTTER 

69.5 

1 

1 

1 

5.5 1 
1 

FLUTTER 

1 

1 

1 

1 

0.0 1 
1 

FLUTTER 

89.3 

1 

1 

6.0 1 

FLUTTER 

1 

1 

0.0 1 

FLUTTER 




the program predicted flutter with a B j[_p of only 10.0 degrees while the 
experiment showed no flutter at all, whereas at 55.1 ft/ sec the program 
predicted no flutter with a of 7.0 degrees while the experiment 
showed flutter with a 6j_g of only 5.5 degrees. Then, for the final two 
velocities the program predicts flutter without any tapping while. In 
the experiment, the wing had to be tapped to 5.5 degrees and 6.0 
degrees, respectively. It should be noted that the program predicted 
the same limit cycle for all B^^.p that produced flutter. 

Note, however, that there la a difference between the two ways 
of "tapping” the wing. In the experiment, the wing was actually tapped, 
and the maximum angle to which this tap would cause the wing to 
oscillate, B^_g, was recorded. In the program, however, the wing was 
Just "locked out" to an Incidence angle, B£_p, and then released. 

Thus, the above comparisons are not entirely valid since the 
aerodynamic history of the experimental wing, before reaching Bi-e> 
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unsteady, whereas the aerodynamic history of the computer wing, before 
reaching was steady. This would certainly have an effect, at 

least for a short time after the maximum angle of attack was reached, on 
the unsteady forces acting on the wing. This effect could, in some 
cases, be great enough to cause a discrepancy between predicted and 
experimental stall flutter onset. 


Reasons for Inaccurate Predictions 


The explanations for experimental and predicted flutter 
discrepancies to follow are qualitative, since the complicated nature of 
this problem makes a quantitative analysis virtually impossible. It is 
hoped that this discussion will aid the reader in understanding both the 
stall flutter phenomenon and the problems Inherent in dynamic stall 
testing and prediction. 


Experimental Measurement Error 

A close look at experimental values and their error estimates 
reveals, immediately, that the discrepancies between them and predicted 
values is not due to experimental measurement error*. For Instance, 
measurement error does not explain the fact that predicted data for 
oscillations about the 1/4 chord elastic axis are consistently low 
whereas predicted data for those about the 37. SZ chord elastic axis are 
consistently high. Bather, <xie would expect a more "normally" 
distributed error. In addition, it gives no hint as to why the 


See Appendix B for error estimation. 


amplitude slopes for the 1/4 chord elastic axis are never correctly 
predicted* 

Inaccuracy In the Wing and 
Wing-Support Construction 

Due to the less-than-perfect fitting of the wing Into the 
frame's wing supports, the wing was able to oscillate vertically with 
amplitudes of up to approximately 1*0 mm. MIT found that, when an 
airfoil was allowed both translational and torsional freedom during 
stall flutter, the airfoil oscillated vertically with anplltudes varying 
between 0.1 mm and 1*0 mm and that the torsional amplitude Increased by 
roughly 2.0 degrees and showed a positive Increase with respect to V^. 
Therefore, It Is entirely possible that this seemingly Insignificant 
free-play, which was strictly a mechanical error In the wing support, 
could have caused, at least In part, the differences between predicted 
and experimental amplitude data* However, this Inaccuracy could not 
have been the only factor Involved In causing these differences since 
MIT reports a positive dui/dV^ value for this two-degree-of-freedom 

motion, whereas PSU finds a dui/dV value equal to zero* Thus, this 

00 

explanation remains only a possibility, at best, until further data Is 
collected with this second degree-of -freedom completely eliminated. 

Flow Interference due to the Frame 
and Wing Supports 

Free-stream flow Interference caused by the proximity of the 
wing-support frame to the wing itself are estimated to be negligible* 
Tuft flow visualization showed that the large endplates eliminated any 
gross flow disturbances that may have been created by the frame's metal 
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bars. 

The aluminum support bars and accelerometer attached to, and 
oscillating with, the wing may have changed the wing's static moment and 
drag characteristics enough to cause significant error. Since these 
attachments do not stall In the same manner as a NACA 0012 airfoil, the 
UTRC program could not properly correct for their preeence. A better 
way to take this effect Into account would be to substitute the static 
data of the wing, with endplates only. Into the UTRC program and then to 
correct for support bar and accel.erometer effects by adding these 
attachments static Cjj and Cj^ values to the wing's unsteady Cjj and Cjj 
values. In this way, the non-linear normal force and moment coefficient 
values peculiar to a dynamically stalling airfoil are calculated for the 

wing only, and any extra forces or moments due to other effects are 
merely added on. 

Inaccurate Static Data 

The recalculated static characteristics, mentioned In Chspter 
IV, did yield, for two recalculated cases, flutter predictions that were 
different from what has been reported herein. However, these 
differences do not significantly change original comparisons with 
experimental results. Furthermore, It Is believed that the greater 
source of error Is due to the manner In which the static data was used 
in the UTRC program, which Is discussed In the previous paragraph, 
rather than to Inaccurate data collection. 
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Tunnel Effects 

Tunnel Interference effects, present during the collection of 
the UTRC data base, could have altered the relationship between static 
and dynamic data especially since the UTRC two-dimensional tunnel has 
such a small test section height compared to the airfoil chord length. 
However, Just how much of an effect this may have had is impossible to 
estimate since the difference between static and dynamic tunnel 
interference effects is currently unkno%m (4, 30). 

Reynold's Ntnber Effects 

The low experimental Reynold's Numbers used in these tests of, 
1.0x10^, to 2.5x10^ could have caused greater oscillation amplitudes 
about the 1/4 chord due to delayed boundary layer reattachment after 
stall. As stated in Chapter II, dynamic reattchment at a Reynold's 
Number of approximately 1.0x10^ occured later than at a Reynold's Number 
of 2.5x10® and caused greater moment hysteresis and negative aerodynamic 
damping. It is, then, very possible that the low experimental values 
used in the PSU experiment, as compared to those used for the UTRC data 
base (9.6x10^), could partly account for the experimental amplitude 
values being higher than those predicted by the UTRC program. However, 
this does not explain why predicted amplitudes underestimate 
experimental amplitudes obtained from the 1/4 chord elastic axis and 
overestimate those obtained from the 37. SX chord elastic axis. 
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Mach WuabT Eff«cti 

The previously unaccounted for discrepancy between 1/4 chord and 
37*SX chord elastic axis data comparisons can be explained In terns of 
dynamic stall type and Mach Nuuber effects* The UTRC data base was 
taken at a free-stream Mach Number equal to 0.325. This relatively high 
subsonic Mach Number would easily cause early, shock-induced flow 
separation and, therefore, prevent the airfoil from entering Into deep 
dynamic stall. This would significantly alter the clockwise moment 
hysteresis loops for the oscillating airfoil, especially at large 
amplitudes where the force and moment overshoots can be large and where 
subsequent dynamic stall behavior Is strongly dependent upon the 
powerful, Fned vortex. In other words, the UTRC program may be 
predicting and 0^^ values characteristic of light stall, when. In 
fact, at lower Mach Numbers, conditions «K>uld have yielded and 
values characteristic of deep stall. This Is borne out also In the 
comparisons done In references 31 and 32. Generally, light stall will 
cause greater negative aerodynamic damping and, therefore, greater limit 
cycle amplitudes. Thus, the UTRC program would be expected to 
overpredict deep stall amplitudes. This Is the case for oscillation 
amplitudes about the 37. SZ chord. One would not expect this effect to 
be quite as pronounced for the 1/4 chord where light stall conditions 
prevail throughout. 
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Three-Dimensional Effects 

rinilly. It mu«t be pointed out thet large errors may result 
from the three-dimensional flow pattern o.er the wing. Figure 30 shows, 
roughly, how the wing stalled statically as visualized with tufts placed 
at the 12.51, 25.01, 37,51, 50.01 and 62.51 chords. The stall reglona 
are Inside of the thin. bUck lines. It can be seen that the wing 
stalls inboard first and that this stall gradually moves outboard along 
the span. This Is Indicative of varying downwash along the span, 
characteristic of three-dimensional wing stall. 

In this study, the three-dimensional static characteristics were 
substituted into a program designed for a two-dimensional alrfoU. It 
is quite possible that this program cannot accommodate such a 
substitution. Perhaps better results would be obtained if the wing's 
two-dimensional static characteristics were used in the program, and the 

three-dimensional effects were taken into account by using a strip 
analysis of the wing. 
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CH/PTER VI 


CONCLUSIONS AND RECOMMENDATIONS 
Conclusions 

In light of the comparisons between predicted and experimental 
stall flutter data and taking the dynamic stall characteristics Into 
consideration, the following conclusions are drawn: 

1. The empirical prediction method developed at UTRC yields 

reasonably accurate values for unsteady, non-llnear C and C 

n in 

for a NACA 0012 airfoil undergoing dynamic stall. 

2. This program, coupled with an accurate dynamics subroutine, 
produces fair estimates of stall flutter amplitude and frequency 
variation with respect to free-stream velocity for a NACA 0012 
wing section undergoing torsional stall flutter. 

3. The large percentage differences between predicted and exp- 
erimental data are probably due to differences between test 
conditions of the PSU experiment and those of the data base 
rather than to program Inaccuracies or random, experimental, 
measurement error. These differences Include Reynolds and Mach 
Number effects, tunnel Interference, vertical oscillations of 
the wing, wing-support bar and accelerometer interference 
effects and three-dimensional flow effects. 
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4« The best empirical prediction method for a given airfoil is 
one which is derived from a data base most closely representing 
the dynamic stall type and unsteady conditions that the airfoil 
in question will experience. 

Recommendations 

Recommendations for further study are as follows: 

1. Reconstruct the wing support so that no vertical free-play 
is possible. 

2. Increase test Reynolds Numbers to correspond to those of the 
UTRC Program data base, 9.6x10^, either by increasing the chord 
length or by using a trip-wire to Increase turbulence over the 
wing surface*. This may decrease C^j vs. cx hysteresis and reduce 
flutter amplitudes about the 1/4 chord. 

3. Attempt to account for low Reynolds Number effects by 
including static hysteresis data in the program data input. 

4. Eliminate all three-dimensional effects by doing the static 
tests in a two-dimensional tunnel and then accounting for three- 
dimensional effects by applying the computer code in a strip 
analysis along the wing span. 


*According to reference 24, this should not have a gross effect on the 
type of dynamic boundary layer separation that the NACA 0012 
experiences. 
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5. Account for effects on flutter predictions due to the 
aluminum, wing-support bars and accelerometers in the following 
way: substitute the static Cj^ and Cj^ characteristics for the 

wing, with endplates only, into the UTRC Program, and then add 
the static C^j and values of the aforementioned attachments to 
the unsteady normal force and moment coefficients of the wing. 


APPENDIX A 


TORSIONAL SPRING CONSTANT, INERTIA 
AND DAMPING RATIO MEASUREMENT 

Torsional Spring Constant 

Torsional springs used in this experiment were constructed using 
a varying number of translational springs connected to the wing via two 
aluminum support bars, as described in Chapter IV. All were linear over 
the range of experimental oscillation amplitudes. Thus, the overall 
torsional spring constant w'ls calculated by measuring the individual 
translational spring const '.nts and by multiplying those values by the 
square of the moment arm. That is, the torque on the wing, T, is given 
by the equation (Figure A-1); 

T»FxR 

Expressing this in terms of the translational spring constant, K.^, and 
assuming that each spring stretches by a value equal to that subtended 
by the circular arc, s, yields the following equation: 

2 2 

T=K^F 0 

This approximation yields only a 1.0 percent error for an oscillation 
amplitude, Aa , of 20.0 degrees. Then, of course, K was calculated 
using the equation 
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Inertia 


Wing inertial values, which varied for different elastic axes, 
were calculated using the relationship 

W'h 

K'W. .a. to b. „ual to bac.u.a of tha vary low .y.t«. 

daapfng, c („. aaotlon). fbo., j ealcul.ta.1 by .ubatftutlo* 

< 04 . Which va. aaaaoad with tha aid of tha oacllloacopa both b.fora and 
aftat aach aaparlnantal run, Into tha abova a,„atlon. Ihaaa valua. ware 
than avaragad ovar all mn. for aach configuration and It wa. tha.a 
final avaragaa that warn uaad In tha UIRC program Input. 

Damping Ratio 

Bafore and aftar aach axparlmantal run. tha wing was "lockad 

out- to an angl. approalmataly tan dagraa. abova It, agulllbrlu. 

position. It wss than ralaasad. and tha attanuatlon curva was 

constructad from tha oscllloscopa (Flgura A-2) by racordlng pask 

acc.lar.tlon valua. a. tha wing urdarwcnt Its dampad oscillations. 

Thasa paaks wara than nonsalUad and curva fitted, a. In Flgura A-3, for 

each different wing Inertia used In the aspcrlnent, according to the 
equation, 
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This was done, with PSU's statistical program, MINITAB, using a linear 
least-squares method. The final value of ^ was calculated using the 
equation 

n 



APPENDIX B 


ACCELEROMETER CALIBRATION AND MEASUREMENT 
ERROR ESTIMATE 

Accelerometer Calibration 

The Endevco 2233 Accelerometer was calibrated as a function of 
torsional oscillation amplitude and frequency. The acceleration of a 
point P (Figure B-1) undergoing sinusoidal, constant amplitude 
oscillations at a distance, R, from the center of rotation Is given by 
the following equation: 

s =-RAoKjo^sin(jt . 

P 

It Immediately follows that the maximum absolute value of the 
acceleration at this point Is 

I s I =RAaa)^ 

I pi 

and, therefore, that the amplitude of this oscillation Is 

Aa*|Sp 1 /Rw^ . 

The next logical step would be to calibrate the accelerometer 
such that its voltage output could be directly correlated to an 
acceleration, in ft/sec^, and to apply this information to the above 
equation. This value, coupled with the knotm radius, R, and frequency 
of oscillation, u , would make It possible to calculate the angular 
amplitude. However, another method was chosen, due to difficulties 
encountered In accurately reproducing low frequency, sinusoidal 
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vibrations, necessary to making this corrrelation. Sixty six data 
points, which gave accurate and consistent visual amplitude and 
“electrical" acceleration values (from the accelerometer output), were 
graphed as in Figure B-2. The vertical axis represents oscillation 
amplitude and the horizontal axis represents peak-to-peak linear 
acceleration, expressed in volte, divided by the square of the frequency 
of oscillation, expressed in cycles per second. These points were curve 
fitted using a linear least-squares , regression line*. Thus, to find 
oscillation amplitudes for any given data point, the accelerometer 
amplitude output was divided by the square of its frequency and matched, 
by the calibration line, to a particular angular amplitude. 

Error Estimation 

Certain amplitude measurement problems were encountered during 
data collection. During the course of the experiment, certain test 
configurations showed highly erratic accelerometer readings* These 
occured either when the flutter amplitude and frequency were so low that 
random low frequency noise was of the same order of magnitude as the 
acceleration or when these dynamic variables were so high that they, 
along with random noise input, overloaded the amplifier. In cases like 
these, visual observation was the only possible means of calculating 
amplitude, and no means of accurately calculating frequency was 
available. However, these incidents were the exception and not the 

*The linearity assumption Is valid since the frequencies only ranged 
from zero to six hertz. 
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rule. Most of the conflguratllone tested yielded fairly good 
correlation between visual and accelerometer measurements. The 
calibration curve of Figure B-*2 la made up of points for which this 
correlation was found to be exceptional. The scatter la most likely due 
to Inaccuracies in reading both the amplitudes from the angular markings 
on the frame and the accelerometer's signals from the oscilloscope. 

Statistically, however, this scatter Is of little consequence. 

It has been calculated, using a PSU statistical program, MINITAB (33), 
that the correlation coefficient between the two graphed variables Is 
0.981. That Is, 98. IZ of the variation In 2M about Its mean, predicted 
from values of 2s/u , Is "explained" by the regression line. A 
correlation coefficient of 1.0 Implies that all the points land right on 
the regression line. In addition, a 99. 9Z confidence Interval for the 

line's slope, which eqtials 15.3 deg'%yc^/sec^, as calculated using a t~ 

2 2 

confidence table. Is ^.64 deg-cyc /sec , assuming that the regression 
line runs through the origin. This means that one can be 99. 9Z 
confident that the slope of the line lies between the Interval shown In 
Flguro B-2. That translates to an error of only 4.2Z. 

This relatively small error was not, then, the major source of 
experimental uncertainty. The largest error contribution csme either 
from the cases, where accelerometer readings were too erratic to record 
or where visual and accelerometer data did not fall to within 2.0 
degrees of eachother (a very rare occurence). This two degree limit was 

*Thls calculation is valid since the number of data points used, 66, Is 
approximately 13Z of the whole and 10. OZ Is required for validity. 
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estimated, rather arbltrarly by the author to be the maximum possible 
error due to Inaccurate ’’eyeballlng'* of oscillation amplitudes* It Is, 
therefore, believed. In view of the above discussion, that this "visual** 
error ^s the limiting factor for experimental accuracy, and, for this 
reason, a nominal error of ^2.0 degrees Is chosen as an estimate of 
amplitude measurement scatter. 

Since frequency could be read directly off the oscilloscope, 
with no subsequent calibration required, error Is limited to the 
accuracy and care with which this task was isost often carried out. In a 
typical case, where the sine wave peaks were roughly three divisions 
apart (on the scope screen) and where each peak could easily be read to 
within a tenth of a division of Its true position, the percent error Is 
calculated to be 6.7 percent. It Is felt that this value Is probably a 
slight overestlmatlon of the typical error and that an overall value 
estimate of approximately 5.0Z Is acceptable. 
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